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I. INTRODUCTION 



There exists by now considerable experimental evidence that Quantum Chromodynamics 
(QCD) is the theory of strong interactions. 
The QCD Lagrangean 

C = -\g^ v + E i>f - m f ) * t (1) 
4 / 

has the simplicity and the beauty of a fundamental theory. 

The notation in ([[]) is standard. QCD is a gauge theory with gauge group SU(3), 
coupled to the quark fields ipf which belong to the fundamental representation {3}. / 
indicates flavour species (u, d, c, s, t, b). 

G^v is the field strength tensor, belongs to the adjoint representation {8}, and has 8 
colour components G a v . 

In terms of the gauge field A a ^ 

G% = d,Al - d v Al + gt hc A\Al (2) 

jabc are structure constants of the gauge group, g is the coupling constant, and a s = 
g 2 /4<K. G^G^ in (HD is a short notation for £Li G^G^, 

Dfi — — igTAfj, (3) 
is the covariant derivative. T a are the generators in the fundamental representation 



rpa rpb 



if abc T c (4) 



with normalization 



Tr (T a T h ) = ^5 ab (5) 
An alternative notation that we will use is 

A^A,.f=Y, A l Ta (6) 



a 



and 



In this notation 



G^ = G^-f = d^A v - d u A^ - ig [Ap, A v \ (7) 



— -G^uGtiu — — -Tr (G^ V G^ V ) (8) 

Tr is the trace on colour indices. 

Most of the evidence for the validity of QCD is obtained by probing short distances. A 
well known property of the theory is indeed asymptotic freedom : the effective coupling 
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constant vanishes logaritmically at short distances, or at large momentum transfers. The 
rate of the decay it — > 77, which is related to triangle anomaly, tell us that the number of 
colours of the quarks is 3. The structure functions of inclusive lepton - hadron scattering 
in the deep inelastic region are consistent with the scale evolution of the coupling constant 
governed by the perturbative beta funcion [|T]-[3|] 

da a (/x) fa 2 (3i 3 

with 

A> = 11 - \N f (10) 
19 

Pi = 51 - -N f (11) 
Nf is the number of flavours with rrif <C \i. The solution of (§) at one loop is 

c,tf) = °' (A2> (12) 
1 + /?o« s hi — 

At large distances the coupling constant becomes large (infrared slavery), and the pertur- 
bative expansion sensless. 

In addition to that the perturbative series is intrinsecally divergent. One could think of it 
as an asymptotic expansion of some funtion of the coupling constant. Usually the existence 
of such a function is proved by resummation techniques, like Borel resummation. For QCD 
this is problematic Q. 

Alternative methods to perturbation theory must be used to explore large distances. 

A fundamental problem in QCD , which involves large distances is confinement of colour. 
At small distances gluons and quarks are visible as elementary constituents of hadrons. 
However neither a free quark nor a free gluon has ever been observed in nature. Upper 
limits have been established for the production of quarks in high energy collisions. I will 
instead quote a limit deduced in the frame of standard cosmological model ||. Assuming 
that the effective mass of quarks in the primordial quark gluon plasma is of the order of 
a few GeV, the ratios n q /n p and n q jn^ of relic quarks and antiquarks to relic protons and 
photons can be computed. One finds n q /n p ~ 1CT 12 to be compared to the experimental 
upper limit n%jn p < 10~ 27 obtained from Millikan - like experiments @. 

In the language of field theory confinement corresponds to the statement that asymptotic 
states are colour singlets. If QCD is the correct theory of hadrons, then confinement of 
colour must be built in it, i.e., it must be derivable from the lagrangean (|]). The problem is 
then to understand by what mechanism asymptotic states of QCD are forced to be colour 
singlets. 



II. LATTICE QCD 

The most successful non perturbative approach to QCD is the lattice formulation [0. 
The idea is to regularize the theory by discretizing the Euclidean space time to a cubic 
lattice with periodic boundary conditions. 
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A generic field theory is defined in terms of the Feynman path integral 



Z [J]= J 



Dip exp — S[(p] — / Jipd, 



x 



(13) 



Z[J] generates all correlation functions. 

The integral is performed on the field configurations <p(x), which are a continuous infinity: 
it is a functional integral. The way to define a functional integral is to discretize x, and then 
go to the limit in which the density of discrete values goes to infinity. Lattice corresponds 
to a step of this discretizing procedure: on a lattice the number of integration variables is 
finite, the integral (0) becomes an ordinary integral and can be computed numerically, e.g. 
by Montecarlo methods. 

I will not survey in detail the lattice formulation of QCD , but only make a few general 
remarks on it ||. 

i) On a lattice computations are made from first principles: the theory is simulated in 
its full dynamical content. 

ii) To go from lattice to continuum both the ultraviolet cut-off (the lattice spacing a) and 
the infrared one (the lattice size La) must be removed. Asymptotic freedom helps to 
solve the first problem: as the coupling constant approaches zero the lattice spacing 
becomes smaller and smaller in physical units and the coarse structure of the lattice 
less and less important. In the usual notation (3 = 2N c /g 2 (N c = 3 is the number of 
colours) 



Al is the physical scale of the lattice regularization, the analog of the usual Ar^s of 
the MS scheme. 

In addition at large /3's scaling is governed by the perturbative beta function: this 
regime is called asymptotic scaling, and helps in extracting physics from numerical 
simulations. 

To eliminate the influence of the infrared cut-off the lattice size must be larger than 
the correlation length. 

Notice that the result of numerical simulations are regularized amplitudes, which have 
to be properly renormalized to get physics. 

iii) An advantage of the lattice formulation is that the theory is quantized in a gauge 
invariant way. The building block is the link U^(n), which is the parallel transport 
from the site h to n + fi 



The measure of the Feynman integral is the Haar measure on the gauge group d£7 M (n) : 
if the group is compact the integration is finite. 




(14) 




(15) 
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Z[J] = f U dU »( n ) ex P l-PS(U)] (16) 

The action can be written |7j 

S(U) = - £ Re[n^(n)] (17) 

As a -> 

/?S(C/) ^--G^G^a 4 + G(a 6 ) (18) 

IP^n) is the parallel transport around an elementary square of the lattice (plaquette) 
in the plane /1, v. 

Any other choice of the action which differs from ( |17D by terms 0{aP) is equally valid, 
and is expected to produce the same physics in the limit a — > 0. In the language of 
statistical mechanics these actions belong to the same class of universality. 

iv) The equilibrium thermodynamics of QCD at temperature T is obtained by euclidean 
quantization with periodic boundary conditions in time (antiperiodic for fermions), in 
the limit of infinite spacial volume. 

On a lattice of size Ng in the three space directions and Nj< in the time direction, the 
above conditions imply Nt *C Ns and 

N T a(P) = i (19) 

In QCD a deconfining phase transition takes place at a temperature T ~ 150 — 
200MeV: above such temperature hadrons melt into a plasma of quarks and gluons, 
and colour is deconfined. 

The existence of such a transition has been established by lattice simulations ||. 
By now no clear experimental evidence exists of quark gluon plasma: however many 
experiments are on the way [[L(J . 



III. CONFINEMENT OF COLOUR 

A quantity which characterizes the long range behaviour of the force acting between a 
quark- ant iquark pair, Q Q, is the Wilson loop, W(C). W(C) is the trace of the parallel 
transport along a closed path C in space time 

W(C) = Trjp expi^^da^j (20) 

P indicates ordering along the path C. 

Any parallel transport from x\ to £2 transforms as a bilocal covariant under a gauge 
transformation Q(x). If 
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r%2 

Pc(x 1 ,x 2 ) = Pexp L^dx" (21) 

P c (x l ,x 2 ) -> n(x l )P e (x ll x 2 )rt(x 2 ) (22) 

For a closed loop x\ = x 2 , and the trace being cyclic W(C) is gauge invariant. Wilson's 
action ([H]) is an example of Wilson loop, with contour C the elementary plaquette. 

If as C we take a rectangle of sides R in some space direction and T in the time direction, 
let us indicate by W(R,T) the corresponding Wilson loop (fig. 1). 

It can be shown that W(R, T) describes the propagation in time of a Q Q static pair at 
distance R and that 

W(R,T) ~ exp[-TV(R)] (23) 

T — >oo 

with V(R) the static potential. 

For a confining potential at large distances 

V(R) = uR (24) 

and 

W(R,T) -expf-aTi?] (25) 

The dependence (^) is known law, and signals confinement, a is the string tension, 

whose empirical value is 

a = -L GeV 2 (26) 
The area law fl25|) is observed in numerical simulations at /3 < /3 C , i.e. below the deconfining 



transition: above it the string tension vanishes [|TT . 

In summary lattice simulations show that colour is confined at low temperature, and 
that quarks and gluons can exist as free particles above some (deconfinement) temperature. 

Since the early days of QCD it was suggested p2}4T4] that the mechanism by which 
colour is confined could be "dual type II superconductivity of the vacuum". Dual means 
that the role of electric and magnetic quantities are interchanged with respect to ordinary 
supeconductors. The idea was inspired by pioneering work on dual models | 15| . In ref. 



|T5[ it was argued that the force between a Q Q pair could be produced by field configura- 
tions analogous to Abrikosov flux tubes in superconductivity. This mechanism makes the 
energy proportional to the length of the tube as in (|24|) , and explains at the same time the 
string like behaviour of the dual amplitude JTB[. After the advent of QCD the suggestion 
came naturally that the flux tubes could be produced by a dual Meissner effect, squeezing 
the chromoelectric field acting between QQ into Abrikosov flux tubes, in the same way as 
happens to magnetic field in ordinary superconductors. 

Lattice simulations have demonstrated the existence of chromoelectric flux tubes between 



QQ pairs [17], supporting the mechanism of confinement by dual superconductivity. A direct 



test of the mechanism and the understanding of its origin are of fundamental interest. 
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IV. BASIC SUPERCONDUCTIVITY[18] 

Superconductivity is a Higgs phenomenon. The effective lagrangean [|T^] 

C = -\F^ U + (D^ (D^) - V(cp) (27) 

with 

nv) = -A*W + ^(fW) a (28) 

and 

Dfj, = — \2eA^ (29) 

couples the electromagnetic field to a charged scalar field, which describes Cooper pairs. 
L is invariant under U(l) global rotations 

v _> e 2ie > (30) 

If /i 2 in (^3) is negative, V(<£>) has the shape of fig. 2. 

if acquires a non vanishing vev $ and the symmetry is broken down to rotations 

by angles multiple of 2n/q = ir/e, q = 2e being the charge of a pair. 

It proves convenient to parametrize the field tp by its modulus tj) and its phase e iqd 

ip = ije iqe (31) 

Under a U(l) transformation by an angle a 

ib -> ih e ^ 6 + a A u ^ A u + d u a (32) 
t/(i) u(i) n n v v j 

The covariant derivative of p becomes 

= e [qe [d^ + iq (0^6 - A^) tfj] (33) 

or, putting 

— Aft — 8^6 (34) 
= e iqe [d, - iglj i> (35) 

A^ is gauge invariant and 

= d„A u - dvAp = F M „ (36) 
In terms of the news fields the effective lagrangean becomes 
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C = -\f^ v + d^dyri) + gVAA - A*V 2 - ^ 4 (37) 
and the equation of motion for the electromagnetic field 

dpFp, + q 2 ^ 2 A u (38) 
Neglecting vacuum fluctuations, ip 2 = ip 2 (its vev) and ( |38l) reads 

+ m 2 i, = (39) 

with 

m = V2qip (40) 



The parametrization fl3"4|) for the field corresponds to the well known fact that, in the Higgs 
phenomenon, the phase of the Higgs field supplies the longitudinal component of the photon, 
when it becomes massive. 

In the gauge A Q = a static configuration, 8qA = 0, doip = implies Ei = F 0i = and 
(|39|) for the space components reads 

diFij + m 2 Aj = 

or 

V AH + m 2 A = (41) 

j = m 2 A is known as London current: physically it implies the existence of a steady current 
at zero E, or, since pj = E, zero resistivity. 

Taking the curl of both sides of ( f4"T|) we obtain 

V 2 H - m 2 H = (42) 

Eq.([42"D means that H penetrates inside the superconductor by a length Ai = 1/m and is 
otherwise expelled from the bulk of it: this fact is known as Meissner effect. 

The key parameter, both for producing zero resistivity and Meissner effect is $ = (\ip\) ^ 
being m 2 = 2g 2 |$| 2 (gg). 

$ 7^ indicates that U(l) symmetry is spontaneously broken. Indeed if the vacuum |0) 
were U(l) invariant, it would be 

U(l)\0) =e iqoa \0) (43) 

and, for any operator O 

(0|C|0) = (0|f/(l) t Cf/(l)|0) (44) 
If O carries charge q, WOU = e iqa O and eq.(|0]) gives 

(0|O|0) = e i9Q (0|C|0) (45) 
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which implies either q = or (0|C|0) = 0. A non zero vev of any charged operator, like ip, 
implies that vacuum is not invariant under £7(1), which is a symmetry of C, i.e. that U(l) 
symmetry is spontaneously broken. 

A non zero vev of any charged operator signals superconductivity: the effective la- 
grangean for that operator, due to the universal coupling to A^, does indeed generate a 
mass for the photon. The ground state of a superconductor is a superposition of states with 
different charges (Cooper pair condensation). 

In the above equations two scales of length appear: A 2 = 1/yU, the inverse mass or 
correlation length of the Higgs field, and Ai = 1/m, the inverse mass of the photon or 
penetration depth of the magnetic field. A superconductor is called type II if Ai ^> A 2 , 
otherwise it is named type I. The formation of Abrikosov flux tubes at intermediate values 
of the external field H, is energetically favorable in type II superconductors. In type I there 
is complete Meissner effect below some value of H c of H, and complete penetration of it 
above H c , with destruction of superconductivity. 



V. MONOPOLES 

A. Generalities on monopoles 

The difficulty in the construction of operators with non zero magnetic charge in terms 
of the gauge fields A^, stems from the fact that monopole configurations have non trivial 
topology. Here we shall briefly review the definition and the classification of monopoles. 



For a more detailed treatment we refer to re. ||20|| , from which most of what we say here is 
extracted. 

The most general form of Maxwell's equation in the presence of both electric (j^) and 
magnetic (if?) current is 



d^ u = j u (46) 
[) F* = i M 



where F* = \e^ upa F pa is the dual of F pv . Free field equations (j u = jf 1 = 0) are trivially 
invariant under the group of transformations 

Ffj, u -> cos OF^ + sin 6F* U (47) 
F^^cos^-sin^ 

or 

E —> E cos 9 + H sin 9 (48) 
H ^ Hcos9 - Esin9 

In particular this holds for 9 = n/2 or 

E -> H H — > —E (49) 

(031) is known as duality transformation. 
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No particles with magnetic charges have ever been observed in nature, in spite of the 
many attempts to detect them. As a consequence Maxwell's equations are usually written 



d,F^ = f (50a) 
d^F*^ = (50b) 



The general solution of eq. (|50b| ) is 

Fnv = d„A v - d v A„ (51) 

with A^ an arbitrary vector field. For any A^ eq.'s ( |50b|) are identically satisfied. They are 
indeed known as Bianchi identities. The very possibility of introducing A^ relies on the 
absence of magnetic charges, eq.'s ( |50b| ). Explicitely they read 

dH 

divH = — - V A £ = (52) 
dt v ; 

If a monopole exists, and we insist in describing the system in terms of A^, the monopole 
must be viewed as one end of a long, infinitely thin solenoid, bringing the magnetic flux to 
infinity, in order to satisfy the first of eq. (|5*2"|) [^TJ(fig. 3). 

The solenoid is known as Dirac string. In order to make this string physically invisible, 
the parallel transport of any charged particle on any closed path around it must be equal to 
1, or 



2?m = e Adx = e$(H) 
J c 



where $ is the flux of the magnetic field across the string. On the other hand by construction 
4vtQm = $ hence 

71 

Qm = Ye (53) 

This is the celebrated Dirac quantization for the magnetic charge. 

In conclusion a monopole field can still be described in terms of A^, provided eq.(|53"D is 
satisfied, at the price of introducing a nontrivial topology. The argument can be generalized 
to non abelian gauge theories and to a generic distribution of charges contained in a finite 
region of space. The idea is to look at the multipole expansion of the field at large distances, 
and to allow for a nonzero magnetic monopole component. We shall restrict for simplicity 
to configurations with zero electric field at some time, say t = 0: a sufficient condition for 
that is that in some gauge 

A {t) = -A (-t), A(t) = A{-t) 

or 

dA 

A = — = 



dt 



at t = 0. Then 
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F 0i = d A l - diA - ig [A Q , A,] = 

Let us denote by A = (A r , Ag, A v ) the components of A in polar coordinates. A r can be 
made zero by a time independent gauge transformation, which does not affect Aq. Let 
A(r, 9, ip) be the parallel transport to infinity along the radius 

A = P exp (i J™ A r (r, 9, <p) dtp) (54) 

A is unitary, A^ = A -1 . Moreover 

<9A 

— + iA r A = (55) 
Or 

Multiplying by A* to the left ( |55D gives 

A' = -iAt— + AU r A = (56) 
or 

A' r is in fact the gauge transformed of A r under A. We will be interested at the behaviour 
of the field at distances r > R and therefore we do not worry about singularities at r = 0. 
We are thus left with A = (0, Ag, A v ), and we look for a configuration behaving as 1/r as 
r — > oo or 

Again we can make ag = by a procedure similar to the one used to make A r = 0. We 
operate a gauge transformation, independent of t and r 

A' = Pexp (i J\g(6 1 ip)d6 ) j (57) 

We are then left with a v alone and the only non vanishing component of F^ u is F dlp = dga^. 
The field equations, outside the space occupied by matter (r > R) read 

d lt JgF>» + [A'>,JgFi»]=a (58) 
g is the determinant of the metric tensor and ^fg = 1/r 2 sin 2 9. For v = tp eq.([58"D gives 

dg^—dga v = (59) 
sin 9 

which has the general solution 

a v = Q(ip)(a + bcos9) 
If we want no singularity at the north pole a v (0) = and 

a v = Q((p)(l-coa6) (60) 
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The equation with v = 9 reads 

d^F* e +[a*,JgF* e 



(61) 



The term with the commutator cancels and the net result is 



d v Q(<p) = (62) 

or Q = cost. The non abelian monopole field is the abelian field times a constant matrix Q: 
in eq. Q59D , (|61|) the term with the commutator which signals the non abelian nature of the 
gauge group has disappeared. 

By our choice the Dirac string lies along the axis 9 = tt. 

If we had chosen the string along the axis 9 = we would obtain 

a v = —Q(l + cos6>) 
The two configuration differ by a gauge transformation 

U = e l2Qv 



If we demand U to be single valued 

exp(i4vrQ) = 1 (63) 

which is the Dirac quantization condition if we keep in mind that in our notation incor- 
porates the coupling constant g. Eq . (|63"D gives 



gQu = y ( 64 ) 

with rrii an integer. 

To summarize a monopole configuration is identified by a constant diagonal matrix Q of 
the algebra, with integer or half integer eigenvalues, up to a a gauge transformation. This 



identification is known as GNO (Goddard, Nuyts, Olive) classification |22| . 

In SU (N) Q is traceless, and is therefore identified by N — 1 half integer eigenvalues: 
they are N — 1 charges corresponding to the residual [/(l)^ -1 gauge invariance, under the 
transformations which leave a matrix Q diagonal. 

For configurations containing many monopoles it can be shown that by use of gauge 
transformations the matrix Q can be constructed, which is the sum of the QiS describing 
the single monopoles. 

What matters is the topology, specifically the homotopy structure of the (N — 1) t/(l)'s 
obtained by GNO construction ||22|| . It can also be shown that the stable monopole con- 
figurations are possible only if the first homotopy group of the gauge group, n 1 (G), is non 
trivial |23]], i.e. if the group G is non simply connected. 

Ui(SU(N)) is always trivial: there exist no stable classical monopole configurations 
unless the symmetry is broken down to a subgroup H with non trivial homotopy, like U(l). 
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B. 't Hooft - Polyakov monopoles 



A well known example of classical monopole configuration is the 't Hooft - Polyakov 
25| monopole. The theory [[2(J is an SU (2) gauge theory coupled to a scalar field in the 



adjoint representation (p: since all fields have integer isospin, the gauge group is SO(3) or 
SU{2)/Z 2 . The lagrangean is 

C = 1 -{D^){D, V ) - -G, U G, V - V(<?) (65) 

where D^ip = d^(p — gA^ A(p is the covariant derivative and 

V^ 2 ) = ^ 2 + -yf (66) 

is the most general potential invariant under the gauge group and compatible with renor- 
malizability. 

If /i 2 < 0, acquires a nonvanishing vev (if), and the symmetry group breaks down to 
U(l), which is the group of isospin rotations around the direction of (<p). An explicit static 
solution of the equations of motion with finite energy can then be constructed. 

With a special choice of the gauge the solution is a "hedgehog" : 

4 = -^*^(r) |> = jF(r) (67) 

(the upper index is an isospin index). 

At large distances A(r) — > 1, F(r) — > 1. By continuity (p must vanish at some point, 
which is identified with the location of the monopole. 

The solution can be gauge rotated to a gauge (unitary gauge) in which <p is oriented 
along the 3-axis <p = (0, 0, $): this gauge transformation is regular everywhere except in the 
points where (p = 0, and coincides with the construction of sect. 1. In fact the configuration 
(P A = 0, A r = 

A e = -a(n ± A u) (68) 
9 

Atp = - \{y ■ n)n — V\o 

with n± = (cosip, sin (p, 0); V = (0,0,1); n = f/r. nj_ A V is 9 independent. Therefore the 
gauge transformation Aq which makes Aq = can be written 



It is easy to check that 



f0 

A e = P exp / gA e d9 = exp i9d(n^_ A in 
Jo 

k\n ■ aAg = cr 3 (69a) 
AlA e Ae + iA^Ae = (69b) 



AlA v A9 + iAj^A fl = a 3 ^ (69c) 



'1 — cos 
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The gauge transformation which makes ip diagonal is called an abelian projection. The 
abelian projection coincides with the GNO construction. The matrix Q of eq.(S4) is in this 
model Q = 2a 3 , corresponding to a monopole of charge 2 Dirac units. 



A gauge invariant field strength can be defined [|4 



F^ = 0- - -0 ■ (D^0 A D v 0) (70) 



whit = <p/\<p\. 

Similarly we can define 



B» = ■ A„ (71) 



is not gauge invariant, since is not covariant under gauge transfromations. In fact 
-> WA,JJ + iWdpU. The identity holds 

= (d^B u - d u B^) - -0 (8^0 A d u 0) (72) 

The two terms in eq. (|72|) are not separately gauge invariant; only their sum is. After abelian 
projection the second term drops and is an abelian gauge field, with vector potential B^ 

F, u = d„B v - d v B^ (73) 

For the solution of the form ( pTD F^ v as defined by eq.([72|) obeys free Maxwell equations, 
except in the point where 0=0 



At large distances (eq. (|67|) ) 



d^ u = 



-» 1 f 

E = H = - — 

g r 6 



F^y for this solution is the field of a pointlike Dirac monopole. Note that in the hedgehog 
gauge of fl67|) only the second term of eq.([72|) contributes to the magnetic field. 

C. The abelian projection in QCD 

In QCD there is no Higgs field. However any operator <&(x) in the adjoint representation 



can act as an effective Higgs field |27|] , and can be used to define monopoles. In what follows 
we shall refer to SU(2) gauge group for the sake of simplicity: the extension to S77(3) is 
trivial. 

In the notation of previous section 

$(x) = $(x)a (74) 
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The abelian projection is by definition a gauge transformation which diagonalizes $(x): like 
in the 't Hooft - Polyakov's monopole configuration analized in sect. 3 ,it is singular when 
= 0. Those zeros will be then world lines of monopoles. 
A gauge invariant field strength 

F tiv = (p- - -0 ■ (D^tp A D v tp) (75) 

can be defined, which in fact is the field of a Dirac monopole in the neighbouring of a zero. 
Again, putting = <f> ■ 

= {d^B u - d v Bj - -0 {d^tp A d v 0) (76) 

In a hedgehog gauge ($ = f) the first term in (|76|) does not contribute to monopole charge: 
the second term carries a net flux of magnetic field 47r/g. After abelian projection instead 
only the first term is different from zero. 

The classical construction of GNO can now be given a quantum mechanical extension. 

Indeed, in each configuration appearing in the Feynman integral, monopoles will be 
located in the zeros of the classical version of $(x). The abelian projection will identify the 
GNO matrix defining the monopole charges: monopole charges will be additive, the field 
strength beeing asymptotically abelian, and the field equations linear. 

Viceversa given a classical configuration with pointlike monopoles a patching of the single 
monopole configurations can be performed by gauge transformations, in such a way that the 
monopole charges are additive: the basic fact is that monopole charge is a topological 
property, which can be described by homotopy, and a group product in the set of paths 
around the Dirac strings of the different monopoles can be defined, such that the winding 
numbers are additive. As a consequence the GNO matrix for the configuration will be the 
sum of the matrices for the single monopoles. What really matters is the charge and the 
location of the monopoles. Any operator in the adjoint representation which is zero in the 
location of the monopoles, and is diagonal with Q will identify an abelian projection which 
makes Q diagonal. 

Defining such operator will allow to label monopoles in different classical configurations 
which enter the Feynman integral. 

In conclusion all monopoles are U(l) monopoles. A monopole species is identified by 
an operator in the adjoint representation. Notice that the monopole charge coupled to the 
magnetic field of eq.(|7^) is invariant under the gauge group: such charges can condense in 
the vacuum without breaking the gauge symmetry 

Notice also that if SU(N) gauge symmetry is not broken, monopoles are unstable. 

D. Monopoles and confinement in QCD 

Any operator 0(x) in the adjoint representation defines a monopole species by abelian 
projection: of course different choices for 0(x) will define different monopole species: the 
number and the location of monopoles will indeed depend on the zeros of 0(x). Dual 
superconductivity means condensation of monopoles in the vacuum. 

Popular choices for the abelian projection correspond to the following choices for 0(x) 
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1) 

0(x) = P(x) (77) 

P(x) beeing the Polyakov line 



P(x, xq) = P exp 



ij A (x,t)dt (78) 



P is the parallel transport in time on the closed path from xq to +oo and back from 
— oo to xq. P(x) transforms covariantly under the adjoint representation. On a lattice, 
due to periodic b.c, P(x) is defined as 

P(n, n ) = n U (n, n + iv ) 

i=l 

Nt beeing the number of links in the time direction. 

2) 

0{x) = Fij 

Fij any component of the field strength tensor. 

3) Maximal abelian H|. This projection is defined on the lattice by the condition that 
the quantity 



M = ]TTr {u,(n)a 3 Ul(n)a 3 } 



be maximum with respect to gauge transformation Q(n). In formulae: 

The operator 0(x) in this case is known only in the gauge where M is maximum, but 
non in its covarian form. In that gauge 

0{n) = £ Ul(n)a 3 U,(n) + U,(n - ^U^n - p) 

The question relevant to physics is: what monopoles do condense in the QCD vacuum to 
produce dual superconductivity, if any. 

As discussed in the introduction this question is better adressed on the lattice. From 
the point of view of physics a possibility is that all monopoles, defined by different abelian 
projections, do condense in the vacuum [p7[| . Is this the case? 



An answer to the above question is possible if we develop a tool to directly detect dual 
supeconductivity. Understanding the problem in U(l) gauge theory will be sufficient, since 
the problem reduces in any case to U{\) after abelian projection. 



1-16 



VI. DETECTING DUAL SUPERCONDUCTIVITY 



To detect dual superconductivity one can use a phenomenological approach, which con- 
sists in the detection of a London current: this approach is discussed in the lectures of D. 



Haymaker, in this course [29 



An alternative method, which will be discussed here, is to detect a nonvanishing vev of a 
quantity with nonzero magnetic charge. In the next section we shall construct an operator 
with non zero magnetic charge, which will be used as a probe of dual superconductivity 



A. U(l) gauge group. 



Pure U(l) gauge theory in the continuum formulation is a theory of free photons. On 
the lattice, the building block of the theory beeing the link U^{n) = exp(iev4 M (n)a), and the 
action beeing the plaquette (Wilson's action) interactions exist to all orders in e. Putting 
ft = 1/e 2 a value ft c exists, ft c ~ 1.0114(2) such that for ft > ft c the system is made of free 
photons. For ft < ft c instead the interaction is strong and Wilson loops obey the area law, 
which implies confinement of electric charge. 

In the confined phase monopoles condense in the vacuum, which behaves as a supercon- 
ductor. This has been rigorously proven in ref. |30|], but only for a special form of the action 
(the Villain action), by showing that a magnetically charged operator exists, whose vev is 



different from zero. In ref. |31| the construction has been extended to a generic form of the 
action. 

The basic idea is the well known formula for translation 



e^lx) = \x + a) 
The analog of (|79|) for a gauge field is 



exp 



i / 11(f) • b(x)d 6 x 



1 1(f)) = \A{x) + b{x)) 



(79) 



(80) 



where A(x) play the role of q coordinates, Il(af) = dC/dA(x) play the role of conjugate 
momenta, and |^4(x)) is a state in the Schrodinger representation. In the continuum 



IT(x,t) = F u *(£,t) 



(81) 



If we choose b as the vector potential describing the field of a monopole of charge m/2g, 
with the Dirac string in the direction n 



b(x,y) 



m n A r 
2g r(r — r- ft) 



r = x — y 



then 



fi(y, t) = exp |i J d 3 x II(x, t)b(x, y) 
creates a monopole of charge m Dirac units at the site y and time t. 



(82) 



(83) 
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The magnetic charge operator is 



Q = J d 3 x VH = J d 3 xV • ( V A A) 
The commutator [Q(t), fi(y, t)] can be easily evaluated, using the basic formula 

q + a 



e ipa , q 



giving 



[Q{t),n{y,t)\ = / d 3 xV- (V Ab(x,y))fi(y,t) 



m 



In eq.(^) the Dirac string potential has been subtracted from b(x,y). 
We will compute the vev 

A = (o\KM\Q) 

as a possible probe of spontaneous symmetry breaking of dual U(l), 1.6. cLS cL probe of dual 
superconductivity. 

On the lattice II* = -Imll 01 and the obvious transcription for /i is, after Wick rotation 
to Euclidean space, 



n(n, no) = exp 



lmU - 0i > n o ) # (n' - n) 



(84) 



b l is the discretized version of b, and the factor (3 comes from the 1/ g of monopole charge 
and 1/g of the field. Of course /x depends on the gauge choice for b. A better definition, 
which makes \x independent of the choice of the gauge for b is 



/j(n> ) =exp|-/?^Re [u 0i {n', n )(e fel - 1 



(85) 



which coincides with the previous equation for small values of cV The line of sites corre- 



sponding to the location of the string must be subtracted: n' in eq.(85) runs on all sites 
except the string. 

When computing (fi) whith the Feynman integral 



n[d^]e" s ii(r7, n Q ) 



(86) 



it can be easily shown that a gauge transfromation on b, is reabsorbed by the invariance of 



the Haar measure of integration, if the compactified form (]85 ) is used. 

Eq. ( ISSp coincides with the construction of ref . [|30| when the action has the Villain form. 

A direct measurement of {fi) gives problems with the fluctuations: indeed \i is the expo- 
nent of a quantity which is roughly proportional to the volume, and therefore has fluctuations 
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of the order V 1 ! 2 . It is a well known fact in statistical mechanics that such quantities are not 
gaussian-distributed, and that the width of their fluctuations does not decrease with increas- 
ing statistics: the same problem occurs with the numerical determination of the partition 
function. Therefore, instead of (p) we will measure the quantity 

p(/3) = ^ln<A*> (87) 

which is the analog of the internal energy in the case of partition function. (//) can then be 
reconstructed from eq.(^), by use of the boundary condition (p) = 1 (/3 = 0), obtaining 

(^> = exp p{P'W (88) 
Jo 

If (p) has a shape like in fig.4. 

The expected behaviour of p is depicted in fig. 5. 

we expect a sharp negative peak around (3 C . Notice that (p) is an analitic function of (3 
for a finite lattice: so it cannot be exactly zero above (3 C . Only in the infinite volume limit 
singularities can develop in the partition function and (p) can be zero in the ordered phase 
without being identically zero. 

The result of numerical simulations with Wilson action are shown in fig. 6. 



It is known that the transition at f3 c is weak first order ||32|| : this means that, approaching 
(3 C from below the correlation length increases as in a second order transition, up to a certain 
value (3, j3 c , where this increase stops. For values of (3 < (3 a finite size analysis can be 
performed, (p) in principle depends on (3, on the lattice spacing and on the lattice size 
L, and for a finite lattice is an analytic function of (3. Since the correlation length £ is 
proportional to (f3 c — [3)~ u for (3 < (3 the (3 can be traded with £ and v. 

Finite size scaling occours when a/C C 1, a/I C 1 are not relevant, and 

(p) = <t>(L^((3-(3 c ),L) (89) 

since ^~ 1 / u ~ ((3 C - (3). 

As (3 approaches (3 C (p) must vanish in the limit L — > oo and 

>= ^ gt > iffli! (90 ) 

L~ x l v p is a universal function of L l l v ((3 c — (3). The quality of this scaling is shown in fig. 7. 

If at infinite volume (p) ~ {(3 C — (3) 6 the index 5 and (3 C can be extracted from that 
universal function. By a best fit we obtain 

v = 0.29 ±0.1 (91a) 
P e = 1.0116 ±0.0004 (91b) 

Monopoles do condense in the confined phase (3 < (3 C . This is consistent with exact 
results existing for the Villain action: as explained in chapter 1 the action in the lattice 
is determined by requiring that it coincides with the continuum action in the limit of zero 
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lattice spacing a — > 0. This gives a great arbitrariness in terms of higher order in a: in 
theories like QCD which have a fixed point where a — > 0, those terms are expected to be 
unimportant. In the language of statistical mechanics they are indeed called "irrelevant", 
and models which differ by them are said to belong to the same class of universality The 
phase transition in U(l) is known to be first order, and therefore strictly speaking these 
arguments do not apply. Anyhow our procedure coincides with that of ref. |3(J and our 
result is consistent with them. 

The virtue of the Villain action is that it allows to perform the transformation to dual 
variables and to get a lower bound for (fi) for (3 < (3 C . With the generic action such a 
procedure is not known: our numerical method supplies to this incovenience. 

A duality transformation can also be performed in supersymmetric QCD with N = 2, 
allowing to demonstrate explicitely the condensation of monopoles ||33|| . 

In what follows we shall look for condensation of {7(1 ) monopoles defined by some abelian 
projection: of course there we do not know the form of the effective lagrangean and we only 
can rely on the numerical methods. 



B. Monopole condensation in SU(2), SU(3) gauge theories 



We shall apply the procedure developed in sect . [VI A| for {7(1) gauge theory to SU(2) and 
SU(3): the extension is kind of trivial, since it consists in repeating the construction of the 
{7(1) theory to the {7(l)'s resulting fron the abelian projection. 

We shall restrict our discussion to the abelian projection of the Polyakov loop: in that 
case the abelian electric field strength F 0i defined by eq. (|70|) only consists of the first term. 
If ip is the Polyakov line, its parallel transport in the time direction D ip = and 

Foi = <p- G 0l (92) 

The commutation relation between the field strength operators 

t), F* k (y, t)] = -i5 ab (<%<9 fc - 5 lk d 3 ) <f< 3 >(£ - y) + if abc (A£(x)<% - AfflS*,) 5^(x - y) 

(93) 

are gauge covariant, and in particular the term proportional to 6 ab , i.e. the commutator for 
a = b, is gauge invariant, and comes from the abelian part of the field strength. This implies 
that the operator constructed in analogy with the {7(1) operator using ( |79|) has magnetic 
charge m. For SU(2) gauge theory in the abelian projection which diagonalizes the Polyakov 
line numerical simulation around the deconfinement phase transition have been performed 



34j| . As explained in sect. II this is done by putting the theory on a asymmetric lattice with 
A; < A',. 

Fig. 8 shows the simulation on a 12 3 x 4 lattices: as in the {7(1) case p shows a very sharp 
negative peak at the deconfining phase transition. 

The position of the peak agrees with the known value of the deconfining transition 



[i~0| The conclusion is that the {7(1) symmetry related to monopole charge conservation is 



spontaneously broken, and hence QCD vacuum is indeed a dual superconductor. 

Similar results for the two kinds of monopoles {7(1) x {7(1) in the case of SU(3) are 
shown in fig. 9 
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VII. OUTLOOK AND FUTURE PERSPECTIVES. 



We have established that QCD vacuum is a dual superconductor. More studies are 
needed to understand what abelian projections define monopole condensing in the vacuum, 
and to test 't Hooft guess that all abelian projections are equivalent ]27 . 



Our method of analysis can also be used to study questions like the order of the decon- 
fining transition and the possible critical exponents. 

In many other numerical investigations |36| , |3T| the number density of monopoles has been 
counted across the deconfining phase transition, as well as their contribution to the observed 
quantities: the indication is that monopoles do indeed determine the dynamics (monopole 
dominance). This information is complementary to our result that QCD vacuum is a 
superconductor: it shows that the degrees of freedom involved in this phenomenon play a 
dominant role in dynamics. 

Discussions with my collaborators, especially with G. Paffuti, and with collegues T. 
Suzuki, M. Polikarpov, D. Haymaker, P. Cea, L. Cosmai are aknowledged. 
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Figure Captions 

Fig.l Wilson Loop 

Fig. 2 Spontaneous breaking of U(l) in the Higgs Model. 
Fig. 3 Dirac Monopole. 

Fig. 4 Expected behaviour of (//) in compact U(l) in 4-d. 
Fig. 5 Expected behaviour of p in compact U(l) in 4-d. 
Fig. 6 p for compact U(l) on a 8 3 x 16 lattice. 

Fig. 7 Test of the scaling eq.(90). Data come from lattices of different sizes. 
Fig.8 SU(2): p vs f3 on a lattice 12 3 x 4 
Fig.9 SU(3): p vs (3 on a lattice 12 3 x 4. 
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